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Abstract 

The vortex line density of turbulent superfluid 3 He-B at very low temperature is deduced by 
detecting the shadow of ballistic quasiparticles which are Andreev reflected by quantized vortices. 
Until now the measured total shadow has been interpreted as the sum of shadows arising from 
interactions of a single quasiparticle with a single vortex. By integrating numerically the quasi- 
classical Hamiltonian equations of motion of ballistic quasiparticles in the presence of nontrivial 
but relatively simple vortex systems (such as vortex-vortex and vortex-antivortex pairs and small 
clusters of vortices) we show that partial screening can take place, and the total shadow is not 
necessarily the sum of the shadows. We have also found that it is possible that, upon impinging on 
complex vortex configurations, quasiparticles experience multiple reflections, which can be classical, 
Andreev, or both. 

PACS numbers: 

67.40.Vs Quantum fluids: vortices and turbulence, 
67. 30. em Excitations in He3 
67.30.hb Hydrodynamics in Hc3 
67. 30. he Vortices in He3 
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I. INTRODUCTION 

The context of this work is quantum turbulence [l| at temperatures T <C T c , where T c 
is the critical temperature. In this regime, the viscous normal fluid component and the 
mutual friction can be neglected, and quantum turbulence takes its purest form: a tangle of 
quantized vortex filaments which move in a fluid without viscosity. 

Experiments at these very low temperatures have produced intriguing results in both 3 He 
and 4 He. In 4 He, McClintock and collaborators discovered that quantum turbulence, initially 
generated by a moving grid, quickly decays, despite the absence of viscous dissipation |2j. 
In 3 He-B, Fisher and collaborators found that quantum turbulence, initially confined in a 
small region, spreads in space and decays 0, J]. These and other results raise challenging 
questions to low temperature physicists and fluid mechanicists alike. 

In the case of homogeneous quantum turbulence, the turbulence's intensity is character- 
ized, at least in the first approximation, by the vortex line density L (vortex length per unit 
volume), a quantity which can be measured using techniques such as second sound and ion 
trapping. From the vortex line density, the typical distance between vortices, £ ~ L -1 / 2 , can 
be inferred. The current understanding of quantum turbulence js| at very low temperatures 
is that, at length scales much larger than £, the nonlinear interaction between the vortex 
lines results in partial alignment and polarization, such that, for k <C l/£, the superfluid 
supports an energy cascade from large scales to small scales, which manifests itself in the 
classical Kolmogorov energy spectrum Ek ~ k~ 5 ^ 3 where k is the wavenumber. Numerical 
simulations performed using the vortex filament model 6j and the nonlinear Schroedinger 
equation model [?], S] confirmed the existence of such spectrum. The energy cascade im- 
plies the existence of an energy sink, and the natural question arises as what should be 
this energy sink in the absence of viscous dissipation. The likely energy sink is acoustic: it 
is thought that kinetic energy decreases due to the emission of phonons by Kelvin waves 
0, [h]] . Kelvin waves are helical displacements of vortex filaments which rotate with angular 
frequency uo ~ k 2 . To efficiently radiate sound, u, hence k, must be very large: at the 
length scale of vortex separation, £, sound radiation is negligible. To bridge this gap we 
have to appeal to the existence of a Kelvin wave cascade process which generates smaller 
scales, and, in analogy to the classical Kolmogorov cascade, shifts the energy to the required 
high wavenumbers k. Numerical simulations revealed that vortex reconnections decrease the 



2 



kinetic energy directly [ll|, ll2| and trigger the Kelvin wave cascade |1 



The details of this scenario still need to be properly understood. First of all, the possibility 
has been raised that there is an energy bottleneck between the Kolmogorov cascade at 
k -C l/£ and the Kelvin wave cascade at k ^> 



l/£ [l4, • Secondly, recently experiments 
ly, ll7( suggest the existence of a new form of turbulence: a less structured, one-scale, 
"ultraquantum" turbulence state (also called "Vinen" turbulence [18]), which decays as 
L ~ in contrast to the more structured, multi-scale "semi-classical" quantum turbulence 
which decays as L ~ t~ 3 / 2 (consistently with the k~ 5 ^ 3 energy spectrum). Thirdly, the nature 
of the spectrum of L is still unclear: if we naively interpret L as a measure of vorticity, 
;he spectrum of L should increase with k if E k ~ /c~ 5//3 , but experiments show otherwise 
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2l|. 



Homogeneous turbulence is clearly the most important turbulence problem, but, as men- 
tioned before, there are also experiments in which turbulence is confined in a fraction of the 
experimental cell, that is to say it is inhomogeneous and it can spread in space. Examples 
of inhomogeneous or anisotropic turbulence are turbulence generated by a vibrating wire 
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231 ]. or grid [2J], fork 25], counterflow 26| and rotating counterflow 27|, [28] in 4 He, and 

30] observed in rotating 



the twisted vortex state accompanied by a moving vortex front 29 
3 He-B. Inhomogeneous turbulence may seem less generic than homogeneous turbulence, but 
is equally worth of attention. The reason is that at very low temperatures, in a pure super- 
fluid, the key difference 3J| between classical and quantum fluid behaviour becomes more 
apparent: vortex reconnections are forbidden in a classical inviscid Euler fluid, but can take 
place in a superfluid. 

The experimental study of quantum turbulence would be greatly facilitated if better 
visualization techniques were available. Classical turbulence can be investigated using a large 
variety of methods: ink, smoke, Kalliroscope flakes, hydrogen bubbles, hot wire anemometry, 
laser Doppler anemometry, particle image velocimetry (PIV), etc. On the contrary, there are 
few techniques available in liquid helium; the most used are second sound and ion trapping in 
4 He and NMR in 3 He. A drawback of these techniques is that they only measure quantities 
which are averaged over a large region, and we know from the study of classical turbulence 
that it is important to have local information about fluctuations. Fortunately this problem 
has been recognized: work is in progress to build smaller sensors, and new measurement 
techniques have been developed. In 4 He, at temperatures above 1 K, a major breakthrough 
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has been the implement at ion of the PIV technique using micron-size spheres made of glass 



and polymers 



32[ and solid hydrogen [33j. 



In the more difficult regime of very low temperatures 3 He-B, the Andreev reflection 
technique pioneered at Lancaster has been a major advance in providing experimentalists 
with a tool for studying turbulence. The technique is based on the fact that the dispersion 
curve E = E(p) of quasiparticles is tied to the reference frame of the superfTuid, so, in 
a superfluid moving with velocity v s , the dispersion curve becomes E(p) + p • v s , where 



p is the momentum 34], |35j. Thus a side of a vortex line presents a potential barrier to 



oncoming quasiparticles, which can be reflected back almost exactly becoming quasiholes; 
the other side of the vortex lets the quasiparticles to go through. Quasiholes are reflected or 
transmitted in the opposite way. The vortex thus casts a symmetric shadow for quasiparticles 
at one side and quasiholes at the other, and, by measuring the flux of excitations, one detects 
vortices and infer the vortex line density. A similar problem of interaction of rotons with 
quantized vortices in 4 He and formation of shadows for R + and R~ rotons was considered 
by Samuels and Donnelly 36]. 



A related problem of Andreev reflection within the vortex core was analyzed in Refs. 37, 
(see also the book of Volovik 3^]). The analysis in cited works was concerned with the 



bound states, whereas our concern is the propagation of thermal excitations outside vortex 
cores. 

In a recent paper 40] we have solved analytically the semi-classical equations of motion 
of ballistic quasiparticles in the presence of a single stationary vortex. When extrapolated 
to a disordered vortex tangle, our result is in agreement with simpler order of magnitude 
estimates which have been used 
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35] to infer the vortex line density in turbulence exper- 



iments. The aim of this article is to develop our understanding of the interaction of ballistic 
quasi-particles and vortices by considering more complex vortex configurations. 



II. GOVERNING EQUATIONS 

For the sake of simplicity, we consider the problem of motion of quasiparticles in the 
(x, y)-plane in the presence of N straight vortex lines aligned in the z direction. The kinetic 
energy of a thermal excitation of momentum p measured with respect to the Fermi energy 
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where p = |p|. Hereafter we use numerical values taken at zero bar pressure [41[ for the 
quantities which are necessary to describe the motion of the excitation: the Fermi velocity 
Vp ~ 5.48 x 10 3 cm/s, the Fermi momentum pp = m*vp ~ 8.28 x 10~ 20 g cm/s, the 
Fermi energy ep = pp 2 /(2m*) ~ 2.27 x 10~ 16 erg, and the effective mass m* ~ 3.01m = 
1.51 x 10~ 23 g, where m is the mass of the 3 He atom. 

Let A be the magnitude of the superffuid energy gap. Near the vortex axis, at radial 
distances r smaller than the zero-temperature coherence length £o = hvp/irAo « 0.85 x 
O -5 cm, the energy gap falls to zero and can be approximated by A(r) »s Ao tanh(r/£o) 
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431 ] . Since we are mainly concerned with what happens to the excitation for r 3> 



£o, we neglect the spatial dependence of the energy gap and assume the constant value, 
A = 1.76kpT c « 2.43 x 10~ 19 erg, where hp is Boltzmann's constant and T c the critical 
temperature. 

The intersection of each vortex line with the (x, ?/)-plane is a vortex point. Each vortex 
point moves with the flow field generated by all other vortex points. The i th vortex point, 
located at the position rj(t) = Xj(t)i + generates the following velocity field at the 

point r: 

V ^ r > *) = 9 i Kl u^ -' 1 ^ - + j( x - ( 2 ) 

Z7r|r — Yi{t)\ 

where r = xi + yj, i and j are respectively the unit vectors along the x and y axes, and 
the circulation of the i th vortex is «j = ±k; the + and — signs denote respectively a 
vortex (anticlockwise rotation in the (x, -plane) and an antivortex (clockwise rotation). 
The quantity 

« = — = — = 0.662 x 10- 3 cm 2 /s (3) 
2m m 

is the quantum of circulation in 3 He-B. The velocity field at the point r created by the N 
vortices is 

i=N 

v s (r, t) = ^ Vi (r, t), (4) 



th 



thus the velocity of the i vortex point is 
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In the presence of vortices the energy of the thermal excitation becomes 



E= yV + Ag + p-v^r, t). (6) 

In writing Eq. ([6]), the spatial variation of the order parameter is not taken into account for 
the sake of simplicity. We also assume that the interaction term p • v s varies on a spatial 
scale which is larger than £n, and that the excitation can be considered a compact object 
of momentum p = p(t), position r = r(t), and e nerg y E = E(p, r, t). This gives us the 
opportunity to use the method developed in Ref. {^J, and consider Eq. ([6]) as an effective 
Hamiltonian, for which the equations of motion are 

dr dE(p, r) e p p 



dt dp J t ^ + An m* 



+ v., (7) 



dp dE(p,r) d 

Tt= dr~ = -^ (p - v * } - (8) 

Eq. ([7]) represents the group velocity of the excitation in the velocity field of the vortices. 
Excitations such that e p > are called quasiparticles, and excitations such that e p < are 
called quasiholes. The right-hand-side of Eq. ([8]) is thus the force acting on the excitation. 

Before solving numerically Eqs. ([7j) and ([H]) it is convenient to rewrite them in dimen- 
sionless form. We introduce the following dimensionless variables: 

(9) 
(10) 

V., - ^v„. (11) 

R=(x,y) = (f,f ) = |, (12) 

T = t/t , (13) 

where t = (, pf/ A = 2.9 x 10~ 6 s. The Hamiltonian, Eq. and the equations of motion, 
Eqs. ([7]) and ([8]) then become: 

H(n, R, r) = A a/ (II 2 - l) 2 + A" 2 + — 7r 2 n • V S (R, r), (14) 
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and 



dX , 2 n 2 - 1 m* , 

— = A ; v d U x + — n 2 V sx R, r , 15 

dr ^(n 2 - l) 2 + A- 2 m sxK h y 1 
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dY , 2n 2 - 1) „ m* - , n . 

— = A ; 7 II, + — 7T VL(R, r , 16 

-— ^ n ^^r + n 3/^r , (17) 



m \ dX y dX 
dliy m* 2 / dV sx dV sy 
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ir m ~ V n ^ + n »l#)- < 18 » 
where the dimensionless parameter A is 

A = (19) 

In our numerical calculations we shall assume the value A = 10 3 . Finally, the dimensionless 
superfluid velocity is 

i=N i=N 

271-IR- R.-M 



V.(R, r) = £ V<(R, r) = V |„_ r „ 7 M« [-i(y " + - *i(r))], (20) 



i=l i=l 

where 1^ = 1 for vortices, I\ = — 1 for antivortices, and 

j=N 



= J2 v ^)- (21) 

III. SINGLE VORTEX 

The numerical solution of Eqs. (|T5|) - (|T8|) which govern the trajectory of quasiparticles and 
vortices requires special care due to the absence of dissipation mechanisms. Our preliminary 
investigations revealed that the most commonly used differential equation solvers, such as 
for example the Runge-Kutta fourth order method, are not satisfactory, even using a very 
small time step; in the case of a single vortex, these solvers failed to conserve the total energy 
and the total angular momentum of the quasiparticle by large amounts (10% or more). In 
the case of more complex, time dependent vortex configurations, energy and momentum of 
quasiparticles would not be conserved, but clearly we could not trust our results if the basic 
conservation laws were not satisfied in the simplest case of a single vortex. 

Ideally, to build the conservation law into the numerical scheme, the numerical method 



must be symplectic and conserve phase-space volume [451 ] . Unfortunately the known sym- 
plectic algorithms are geared to problems (mainly in the context of gravity) in which the 
Hamiltonian has the additive form H = T(p) + V(q), where p and q are the generalized 
momenta and positions, T is the kinetic energy, and V the potential energy, whereas in our 



problem the variables p and q appear in nonlinear combinations. The second difficulty is the 
stiffness of our equations of motion, as very rapid time-scales appear at the Andreev turning 
points. After some experimenting, we have found that we can solve the governing equations 
with satisfactory accuracy using the Matlab code odel5s, which is a quasi-constant step 
size implementation of the numerical differentiation formulas (NDF) particularly efficient 
for solving stiff problems (for detailed description of the odelSs Matlab solver and corre- 
sponding software see Ref. [461] ). When solving Eqs. (IT5"|) - (IT8"|) . error tolerances were lowered 
until the particle trajectory had sufficiently converged, in particular at reflections. 

To test our numerical method we determine the trajectories of excitations in the presence 
of a single vortex located at the origin, and compare the results with previous analytical 



results 



40j . The velocity field of the vortex is simply 



V.(R) = ^(-iy+jX). (22) 

Since the vortex does not move, this velocity field is time- independent, and the governing 
equations (I15j) - (jl8j) have two integrals of motion: the first is the energy, defined by the 
Hamiltonian, Eq. ( [141) ; the second is the ^-component of the angular momentum, which is 

J = Tl y X - U X Y. (23) 

The initial conditions at r = for our calculations are the following. The initial mo- 
mentum is IIo = (1.0001, 0) and corresponds to a quasiparticle moving in the x direction 
with energy E = A + ksT for temperature T « 0.1T C . The initial position is (X ,Y ) 
with X = — 10 4 , far away from the vortex. We study the trajectory of the quasiparticle as 
a function of Yq, which plays the role of impact parameter. Fig. [1] shows results for some 
typical values of Y . We distinguish three cases: 



Case 1: For Y > we have no reflection, in agreement with previous work 40] . For 
example, Fig. [2] (left) shows the quasiparticle's trajectory for Y = 10; Fig. [2] (right) shows 
that II 2 — 1 remains positive at all times r, which means that the quasiparticle retains its 
nature of quasiparticle. Fig. [3] confirms that our numerical method conserves energy and 
angular momentum very well. The left hand side of Fig. [3] shows that the relative error in 
the energy, 5h(r) = (H(r)—Ho)/Ho, where Hq is the initial energy, is less than 6 x 10~ 10 ; the 
right hand side of Fig. [3] shows that the relative error in computing the angular momentum, 
Sj(r) = (J(t) — Jo)/ Jo, where J is the initial angular momentum, is less than 2.5 x 10~ 9 . 
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Case 2: If Y < but \Y \ is not too large, the incident quasiparticle is Andreev reflected, 
as shown for example in Fig. H] (left) for Yq = —10. Fig. H] (right) shows that II 2 — 1 
changes sign, thus confirming that, upon reflection, the quasiparticle becomes a quasihole. 
In this calculation, the numerical errors in conserving energy and angular momentum are 
5h < 8 x 10~ 10 and Sj < 2 x 10 -9 respectively. Fig. shows another Andreev reflection, this 
time for Y = -205. 

In our previous paper 40j we determined the distance from the vortex at which, if Yq < 0, 
the incident quasiparticle is Andreev reflected; the dashed-dotted (red) curve in Fig. [1] 
marks this location. It is apparent that there is a maximum value of \Yq\ past which a 
quasiparticle with energy e p is not Andreev reflected; this value (in our dimensionless units) 
is approximately equal to So = 37r(A /e p ) 2 « 269 where we used e p = epijl 2 — 1) ~ 0.0002ep 
for II = 1.0001. We call 5*0 = 269 the (dimensionless) Andreev shadow of a single vortex to 
quasiparticles of that particular (dimensionless) momentum IT. 

Case 3: Finally, if > —S > Y , the quasiparticle's trajectory is deflected by the vortex 
but remains a quasiparticle. 

IV. VORTEX- VORTEX PAIR 

The velocity field of two vortices is time- dependent, thus the Hamiltonian of the thermal 
excitation has no integrals of motion. Unlike the previous case of a single vortex, H and J 
are not conserved. The only quantity in the problem which is constant is the the distance 
between the vortices. There are two cases to consider: two vortices of the same circulation 
(vortex- vortex pair), and two vortices of the opposite circulations (vortex-antivortex pair). 
This section is concerned with the former. 

Two vortices of the same circulation at distance d from each other rotate around a point 
halfway between them with velocity 

= m m 

In dimensionless variables we have 

Vrot = ^=, D=L (25) 
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Far from the vortices, the velocity of the quasiparticle can be estimated from Eqs. (|T5|) 
and (HSJ): 

^ ~ 2A 2 (n 2 - l)n. (26) 

GST 

To get a more clear picture of the problem, it is useful to make the following simple 
estimates. Away from the vortices, the velocity of the quasiparticle is approximately 400. 
If D = 10, the velocity of the vortices is approximately 0.016. For a typical timescale of 
approximately 25 to 30 time units, the vortices travel approximately the distance 0.4 to 0.5, 
which means that they rotate about the center of rotation by an angle 5(f) ~ 2.3° to 3°. If 
D is larger than 10 the vortices move slower and 5(f) is even smaller. We conclude that, in 
the first approximation, the vortex system is static to quasiparticles with the momentum 
n = (1.0001, 0). However, if n = (1.00005, 0), in the corresponding time scale the vortices 
move by a more substantial angle, 5(f) ~ 4.5° to 6°. If D = 10, a significant displacement can 
be observed for quasiparticles with n = (1.00002, 0) and the vortex configuration cannot 
be considered static. 

We proceed with our calculations and consider two vortices of the same circulation at 
distance D = 1000 from each other, initially located at positions (Q x i, Qyi) = (0, —500) 
and (Q X 2, Q y 2) = (0, 500). We integrate the equations of motion of quasiparticles with 
IIo = (1.0001, 0), keeping X® = — 10 4 fixed and varying Yq. In this configuration, the 
relative angle 9 between the direction of the incoming quasiparticle (the x axis) and the line 
which joins the two vortices (the y axis) is 9 = 90°. Fig. [6] illustrates the trajectories of the 
quasiparticle for Y = 316 as well as the path of the vortices; Fig. [7J (left) shows that the 
incident quasiparticle is reflected as a quasihole; since the p • v s term in the Hamiltonian 
is time-dependent, the energy is not constant during the evolution, as confirmed in Fig. [7J 
(right). We find that the Andreev shadow of the first vortex of the pair is Si = 290, slightly 
more than Sq (the Andreev shadow of an isolated stationary vortex), and that the shadow 
of the second vortex of the pair is S 2 = 184, which is slightly less than Sq. Note that in this 
case 5*i + S 2 ~ 25*0. 

If the relative angle 9 between the direction of the incoming quasiparticle and the direction 
between the vortices is reduced from 9 = 90° to 9 = 45°, the Andreev shadow of the first 
vortex decreases from Si = 290 to Si = 272, but the Andreev shadow of the second vortex 
increases from S 2 = 184 to S2 = 205, so that the total shadow S% + S2 of the vortex 
configuration remains approximately equal to twice the shadow S of a single isolated vortex. 
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For example, if the vortices are located at (Q x i, Qyi) = (—353, —353) and (Q x2 , Q y2 ) = 
(353, 353), the angle is 9 = 45°, Si = 272, S 2 = 205, and Si + S 2 = 477. If 9 is further 
reduced from 9 = 45° to 14.1°, the total shadow of both vortices increases and reaches its 
maximum size: Si = 263, S 2 = 242, and Si + S 2 = 505. If 9 > 14.1° the two shadows merge, 
and the two vortices screen each other. If 9 — 0° then S± + S 2 = 251 only. 

Qualitatively, the same behaviour is observed if the distance between the vortices is even 
smaller, e.g. D = 100. If the angle between the direction of motion of the excitation and 
the direction of the line through the vortices is 9 = 90°, the shadow of the first vortex at 
(Qxi, Qyi) = (0, 50) is only Si = 44.7 (significantly smaller than So) and the shadow of the 
second vortex at (Q x2 , Q y2 ) = (0, —50) is S 2 = 443.3 (significantly larger than So), and the 
total shadow is Si + S2 = 488. If 9 is reduced the two shadows merge when 9 ~ 75°; at 
this angle Si + S 2 = 526. At 9 = 45°, the total shadow is Si + S 2 = 510, and at 9 = 0° 
Si + S 2 = 469. 

What happens if D is further reduced? Consider for simplicity the angle 9 = 90°: for large 
D the two vortices have independent shadows (one reduced in size, the other increased in 
size). By decreasing the distance between the vortices the shadows approach each other; for 
distances D < 16 they merge and the vortex configuration has a single shadow independent 
of 9, as if it were a single vortex of strength approximately equal to 2k. 



V. VORTEX- ANTIVORTEX PAIR 



A vortex and an antivortex, set at distance D from each other, move through the fluid 
with (dimensionless) translational velocity 

"—sb- D = h (27) 

The same estimates which we have made at the beginning of the previous section apply. 
As before, we consider quasiparticles with initial momentum n = (1.0001, 0) and initial 
position (X , Y ) with X = — 10 4 fixed and varying Y . Firstly we consider the case in 
which the vortex-antivortex pair and the quasiparticle move in the same direction. Let 
the positive (anticlockwise) vortex be located at (Q x i, Q y i) = (0, 500) and the negative 
(clockwise) vortex be at (Q x2) Q y2 ) = (0, —500), with D = 1000 the separation between 
the vortices. We find that the total shadow of the vortex configuration is Si + S 2 = 774, 
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as shown in Fig. [HI Secondly we consider the case in which the vortex-antivortex pair 
and the quasiparticle move in the opposite directions, letting the positive vortex be at 
(Qxi, Q y i) = (0, —500) and the negative vortex at (Q X 2, Q y 2) = (0, 500): we find that the 
total shadow is greatly reduced: Si + S 2 = 332, as shown in Fig. [9J In both cases, during 
the timescale of the calculation, the vortex pair moves by only 0.01. We conclude that the 
relative motion of vortices and excitations has a strong effect on the shadow. 

Now we examine the dependence of the Andreev reflection on D. Firstly, we consider 
the case in which the quasiparticle and the vortex pair move in the same direction. We said 
that, at D = 1000, the total shadow is S\ + S 2 = 774. If we reduce D, the total shadow 
increases, and at the value D = 940 the two shadows merge into a single shadow of size 
Si + 5*2 = 940. Upon further reduction of D, the total shadow decreases; for example, when 
D = 100, 51 + 52 = 122, and when D = 10 we have Si + S 2 = 34. Secondly, we consider the 
case in which the quasiparticle and the vortex pair move in opposite directions. We have said 
that if D = 1000 the total shadow is only Si + S 2 = 332. If D decreases, Si + S 2 decreases 
too: at D = 100 and D = 10 we have Si + S 2 = 222 and 88 respectively. We conclude that, 
independently of D, the total shadow of a vortex pair travelling in the opposite direction of 
the quasiparticle is about half that of a vortex pair travelling in the same direction. 

The shadow of the vortex-antivortex pair also depends on the angle a between the direc- 
tion of propagation of the quasiparticle and the direction of motion of the vortex-antivortex 
pair. Consider a vortex-antivortex pair of separation D = 100. We have already seen that 
if a = (vortex pair moving in the same direction as the quasiparticle), then S = 122. If 
the positive vortex is at (35.33, 35.33) and the negative vortex is at (—35.33, —35.33), the 
angle is a = —45° and the shadow increases to 168. If the positive vortex is at (50, 0) and 
the negative vortex at (—50, 0), the angle is a = —90° and the shadow increases further to 
S = 233. Finally, as said before, if the positive vortex is at (0, —50) and the negative vortex 
is at (0, 50) (the vortex pair moving in the direction opposite to that of the quasiparticle), 
then a = -180° and S = 222. 

VI. MANY VORTICES 

We assume again that the quasiparticle has initial momentum n = (1.0001, 0) and 
initial position (Xq, Yq) with Xq = — 10 4 ; we vary Yq and determine the total shadow of 
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some simple vortex configurations. 

In the first numerical experiment we initially place five vortex points in the square —250 < 
X < 250, —250 < Y < 250. More precisely, the initial positions of the vortices are (—250, 0), 
(0, 0), (250, 0), (250, -250) and (-250, -250). 

If the vortices have the same (positive) polarity, they rotate around each other, forming 
the 2-dimensional equivalent of a vortex bundle of total circulation 5k; we find (see Fig. [TUj) 
that the total shadow of the vortex configuration is S = 1238, which is less than five times 
the shadow of five individual vortices. 

If we change the sign of some of the vortices, the total shadow which is cast changes 
dramatically. For example, let the vortices at (0, 0), (0, 250) and (250, 0) be positive, and 
the vortices at (—250, 0) and (0, —250) be negative. The net circulation is now k and the 
total shadow is S = 585 (see Fig. [TT1) . which is much less than the previous value, but 
still more than the shadow of a single isolated vortex. A similar value of the total shadow, 
S = 589, is obtained if the positive vortices are at (0, 0), (0, 250), (0, —250) and the 
negative vortices are at (—250, 0), (250, 0), which corresponds to the same total circulation 
k as before, see Fig. [121 

In the second numerical experiment we increase the vortex density and place ten vortices 
(twice as many as before) in the same square —250 < X < 250, —250 < Y < 250. We 
consider three cases: (i) ten vortices of the same polarity (total circulation 10k, see Fig. [T3l) . 
which yields the total shadow S = 2445 (about twice the shadow of a bundle of five vortices); 
(ii) five vortices and five antivortices (total circulation is zero, see Fig. fT4l) . which yields 
S = 902; (iii) five vortex-antivortex pairs (total circulation zero, see Fig. [151), which yields 
S = 209. 

We have found that it is possible that, upon impinging on complex vortex configurations, 
quasiparticles experience multiple reflections, which can be classical, Andreev, or both. An 
example of multiple Andreev reflection is shown in Fig. [T6l 

VII. CONCLUSIONS 

In conclusion, the above numerical experiments with simple vortex configurations show 
that partial screening takes place. The total Andreev shadow of a vortex system is not 
necessarily the sum of the shadows of individual vortices, and depends not only on the 
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distance but also on the relative orientation between quasiparticles and vortex motion. This 
does not mean that the interpretation given to recent experiments is incorrect. It is possible 
that, for a large, random vortex system, the partial screening effects which we have found 
average out. If this is the case, screening effects can be taken into account by introducing a 
prefactor probably of order one for the total shadow, hence for the vortex line density which 
is inferred. Numerical investigations in 3-dimensions with realistic vortex line density are 
needed. 

How random are vortex configurations of current experiments? Probably only the recently 



discovered 
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171 ] ultraquantum regime is truly random. Homogeneous isotropic turbulence 



contains coherent vortex structures 



471 . |48| . |49( and is organized in scales with different energy 



per scale. On the other hand, provided we are interested only in large scale properties 
averaged over a large region, such as the total vortex length, the partial screening effects 
can be accounted as said above. 

3, Q 
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and 



The situation is very different when we move to rotating turbulence 
inhomogeneous turbulence, particularly if there are turbulent fronts. In these cases there is 
large scale anisotropy, and the Andreev reflection technique must be used with more care 
than we used to. The good news is that, by combining Andreev reflection measurements in 
different directions and numerical calculations such as those we have presented, it should be 
possible to gain more information about the geometry and the anisotropy of the turbulence. 

Our results also indicate that the problem of interaction between rotons and quantized 
vortices in 4 He [36|] , leading to calculation of the mutual friction, should be reconsidered in 
view of possible screening effects analyzed above in Sees. [TV] and [V] 
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FIG. 1: (Color online). Trajectories of excitations with initial momentum IIo = (1.0001,0) and 
initial position (Xq,Yo) with Xq = — 10 4 for different values of Yq in the presence of a single 
(positive) vortex at the origin (marked by the dot). The (anticlockwise) superfiuid velocity field 
of the vortex is indicated by arrows. Quasiparticles trajectories are solid (purple) lines, quasiholes 
trajectories are dashed (green) lines. The dashed-dotted (red) curve denotes the locus of Andreev 
reflection. 
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FIG. 2: Left: trajectory of the quasiparticle with initial momentum IIo = (1.0001, 0) and position 
(Xq, Yq) = (— 10 4 , 10) in the presence of a single anticlockwise vortex at the origin. The arrows 
indicate the direction of motion. Right: plot of LT 2 — 1 vs time r corresponding to the left part of 
the figure; note that the quasiparticle remains a quasiparticle. 




FIG. 3: Plot of relative energy variation 5h = {H — Hq)/Hq vs time r (left) and relative angular 
momentum variation 5j = ( J — Jq)/Jq (right) corresponding to Fig. [2j 
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FIG. 4: Left: trajectory of quasiparticle with initial momentum IIo = (1.0001, 0) and position 
(Xq, Yq) = (— 10 4 , —10) in the presence of a single anticlockwise vortex at the origin. The arrows 
indicate the direction of motion. Right: plot of IT 2 — 1 vs time r corresponding to the left part of 
the figure; note that the quasiparticle becomes a quasihole. 
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FIG. 5: Left: trajectory of the quasiparticle with initial momentum IIo = (1-0001, 0) and position 
(Xq, Yq) = (— 10 4 , —205) in the presence of a single anticlockwise vortex at the origin. The arrows 
indicate the direction of motion. Right: plot of LT 2 — 1 vs time r corresponding to the left part of 
the figure; note that the quasiparticle becomes a quasihole. 
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FIG. 6: Left: trajectory of the quasiparticle with initial momentum IIo = (1.0001, 0) and position 
(Xq, Yq) = (— 10 4 , 316) in the presence of a vortex-vortex pair; the first vortex is at (Q x i, Qyi) = 
(0, —500), and the second vortex is at (Q X 2, Qyi) = (0, 500). The solid arrow indicates the 
incident quasiparticle, the dashed arrow the reflected quasihole. Right: Directions of motion of the 
two vortices. 




FIG. 7: Left: plot of LT 2 — 1 vs time r corresponding to the left part of Fig. [6j note that the 
quasiparticle becomes a quasihole. Right: plot of relative energy difference 5h = (H(r) — Hq)/Hq 
vs time r. 
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FIG. 8: Andreev reflection of quasiparticles from the vort ex- ant i vortex pair travelling in the same 
direction. The solid lines denote quasiparticles travelling left to right; the dotted lines denote the 
reflected quasiholes. The solid double arrow denotes the path of the vortex, initially located at 
(Qxi, O y i) = (0, 500) with separation D = 1000; the dash-dotted double arrow denotes the path 
of the antivortex, initially located at (Q X 2, Qy2) = (0> —500). The thick vertical grey lines denote 
the shadows of the vortices. 
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FIG. 9: Andreev reflection of quasiparticles from the vortex-antivortex pair, with separation D = 
1000, travelling in the opposite direction. The solid lines denote quasiparticles travelling left to 
right; the dotted lines denote the reflected quasiholes. The dash-dotted double arrow denotes the 
path of the vortex, initially located at (Q x i, Q y i) = (0, —500); the solid double arrow denotes the 
path of the antivortex, initially located at (Q X 2, Qy2) — (0, 500). The thick vertical grey lines 
denote the shadows of the vortices. Note that the shadows are much smaller than in Fig. 
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FIG. 10: Trajectories of quasiparticles (solid lines), of initial momentum IIo = (1.0001, 0) and 
initial position (— 10 4 , Yq) for varying Yq, interacting with a configuration of five positive vortices 
(denoted by crosses) initially located at (-250, 0), (0, 0), (250, 0), (250, -250) and (-250, -250). 
The reflected quasiholes are indicated as dotted lines; the shadow of the vortex configuration is the 
thick grey line. 
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FIG. 11: Trajectories of quasiparticles (solid lines), of initial momentum IIo = (1-0001, 0) and 
initial position (— 10 4 , Yq) for varying Yq, interacting with an initial configuration of three positive 
vortices (denoted by crosses) located at (0, 0), (250, 0) and (0, 250) and two negative vortices 
located at (—250, 0) and (0, —250) denoted by circles. The reflected quasiholes are indicated as 
dotted lines; the shadow of the vortex configuration is the thick grey line. 
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FIG. 12: Trajectories of quasiparticles (solid lines), of initial momentum IIo = (1.0001, 0) and 
initial position (— 10 4 , Yq) for varying Yq, interacting with an initial configuration of three positive 
vortices (crosses) located at (0, 0), (0, 250) and (0, —250) and two negative vortices (circles) located 
at (250, 0) and (—250, 0). The reflected quasiholes are indicated as dotted lines; the total shadow 
of the vortex configuration, S = 589, is the thick grey line. Note the gap in the shadow within the 
interval -95 < Y < 0. 
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FIG. 13: Schematic trajectories of quasiparticles (solid arrows), of initial momentum IIo = 
(1.0001, 0) and initial position (-10 4 , Y ) for varying Yq interacting with an initial configura- 
tion often positive vortices (crosses) located at (-250, 0), (-250, 250), (-125, -125), (0, -250), 
(0, 250), (125, 125), (250, -250), (250, 0) and (250, 250). The reflected quasiholes are indicated 
as dashed arrows. The thick grey vertical lines indicate the shadows of the vortices. The total 
shadow is S = 2445. 
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FIG. 14: Schematic trajectories of quasiparticles (solid arrows), of initial momentum IIo = 
(1.0001, 0) and initial position (— 10 4 , Yq) for varying Yq, interacting with an initial configuration 
of five positive vortices (crosses) located at (-250, -250), (-250, 250), (-125, -125), (0, 250), 
and five antivortices (circles) at the (-250, 0), (0, -250), (125, 125), (250, -250) and (250, 250). 
The reflected quasiholes are indicated as dashed arrows. The thick grey vertical line indicates the 
shadow of the vortices. The total shadow is S = 902. 
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FIG. 15: Schematic trajectories of quasiparticles (solid arrows), of initial momentum IIo = 
(1.0001, 0) and initial position (— 10 4 , Yq) for varying Yq, interacting with an initial configuration 
of five vortex-antivortex pairs located at: vortex at (—250, —240) and antivortex at (—250, —250), 
vortex at (—250, 250) and antivortex at (—250, —240), vortex at (—125, —30) and antivortex at 
(—132, —36), vortex at (0, 115) and antivortex at (5, 107), vortex at (250, —125) and antivortex 
at (250, —135). (Note that the distance between the vortex and the antivortex in each pair is 
10 non-dimensional units and cannot be distinguished on this figure.) The schematic reflected 
quasiholes are indicated as dashed arrows. The thick grey vertical lines indicate the shadows of 
the vortices. The total shadow is S = 209. 
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FIG. 16: Left: trajectory of the quasiparticle with initial momentum IIo = (1.0001, 0) and position 
(Xo, Yq) = (— 10 4 , 10) in the presence of five positive vortices and five negative vortices as in 
Fig. Q31 Note the multiple reflections before the particle's escape from the vortex region. The 
arrows indicate the direction of motion. Right: plot of LT 2 — 1 vs time r corresponding to a small 
time interval just before the escape. Note that the quasiparticle turns into a quasihole many times, 
before escaping as a quasiparticle. 
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